Optimization of quantum interferometric metrological sensors 
in the presence of photon loss 
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We optimize two-mode, entangled, number states of light in the presence of loss in order to 
maximize the extraction of the available phase information in an interferometer. Our approach 
optimizes over the entire available input Hilbert space with no constraints, other than fixed total 
initial photon number. We optimize to maximize the Fisher information, which is equivalent to 
minimizing the phase uncertainty. We find that in the limit of zero loss the optimal state is the 
so-called NOON state, for small loss, the optimal state gradually deviates from the NOON state, 
and in the limit of large loss the optimal state converges to a generalized two-mode coherent state, 
with a finite total number of photons. The results provide a general protocol for optimizing the 
performance of a quantum optical interferometer in the presence of photon loss, with applications 
to quantum imaging, metrology, sensing, and information processing. 

PACS numbers: 42.50.St, 42.50.Ar, 42.50.Dv, 42.50.-p 



Quantum states of light play an important role in 
applications including metrology, imaging, sensing, and 
quantum information processing [l]]. In quantum inter- 
fcrometry, entangled states of light, such as the maxi- 
mally path-entangled NOON states, replace conventional 
laser light to achieve a sensitivity below the shot-noise 
limit, even reaching the Heisenberg limit, and a resolu- 
tion well below the Rayleigh diffraction limit @- For 
an overview of quantum metrology applications see, for 
example, Ref. [l]. However, for real- world applications, 
diffraction, scattering, and absorption of quantum states 
of light need to be taken into account. Recently it 
has been shown that many quantum-enhanced metrol- 
ogy schemes using NOON states performpoorly when a 
considerable amount of loss is present [!, 0, Hj]. How- 
ever, our team has also discovered a new class of entan- 
gled number states, which arc more resilient to loss Q. 
These so-called M&M' states still outperform classical 
light sources under a moderate 3 dB of loss. 

In this work, we systematize the numerical search for 
optimal quantum states in a two-mode interferometer in 
the presence of loss. We employ the Fisher information 
to obtain the phase sensitivity of the interferometer. An 
exhaustive review and application of the Fisher informa- 
tion concept to the sensitivity of a March- Zchndcr in- 
terferometer, particularly in the zero loss case, has been 
presented in the recent work by Durkin and Dowling 0] . 
The chief utility of the Fisher information approach is 
that it provides a bound on the phase sensitivity, even 
in the absence of a fully specified detection scheme, and 
is now widely adopted in studies of interferometer sensi- 
tivity. Such numerical optimization has been previously 



carried out in the absence of loss, and with loss over a 
restricted class of input states 0, Here, we provide 
a completely general optimization scheme that is applied 
to the two-mode interferometer, but also has application 
to the optimization of linear optical systems for quantum 
linear optical information processing 
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Using this scheme, we first recover the well-known fact 
that NOON states are optimal in the absence of loss [![. 
For large loss, the optimal states belong to a class of two- 
mode coherent states with finite total photon number. 
The optimization procedure yields the optimal Fisher in- 
formation - and hence the minimal phase uncertainty - 
for every level of loss. The validity of our numerical opti- 
mization is verified using several methods, including ge- 
netic algorithms and simulated annealing, and the close 
agreement among these methods provides evidence that 
we are indeed finding the global optimum. 

In quantum optics, photon loss is typically modeled by 
a beam splitter that routes photons out of the interfer- 
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FIG. 1: (Color online) Abstract interferometer condensing 
the input state plus the first beam splitter into the first box, 
followed by two propagating modes with loss modeled by ad- 
ditional beam splitters. The box on the right includes a beam 
splitter and the photon-number resolving detectors. 
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ometer 12[ . In implementing this model, we first enlarge 
the Hilbert space to include modes representing the scat- 
tered photons and then, after the scattering, trace out 
these modes. Here we extend the beam splitter model 
for photon loss to two propagating modes that represent 
the two paths in the optical interferometer (Fig. |T|) . 

For an interferometer with two input ports A and B 
as depicted in Fig. [T] an arbitrary pure-state input with 
N photons can be written as 1^) input = Ylk=o c k\N — 
k, k)A,B where Ck are the input amplitudes to be opti- 
mized. The ket \N — k,k)A,B is a basis state in which 
N — k and k photons are in mode-A (detection arm) 
and mode-B (control arm), respectively. Such a quan- 
tum state resides in an N + 1 dimensional Hilbert space. 
The transformation of the quantum state by any pas- 
sive lossless optical elements, such as beam splitters, 
phase shifters, and mirrors, can then be described by 
an (TV + 1) x (N + 1) unitary matrix. 

When the propagation channels suffer from photon 
loss, we need to consider the total density matrix that 
includes all the scattered photon modes. Then we obtain 
the reduced density matrix for the two interferometer 
modes by tracing out these additional modes [~L3| . This 
reduced density matrix is associated with a much larger 
Hilbert space of dimension (N + l)(N + 2)/2, which in- 
cludes all states with a total of N, N — 1, photons in 
the two interferometer modes. 

For the NOON state, combined with its specific detec- 
tion scheme, the density-matrix description of photon 
loss may be avoided and the state- vector approach can be 
adopted @, 0, U | ■ Previously, we have used the density- 
matrix approach to lossy interferometers for particular 
input states of light (namely, the MhM' states 0]). Our 
method here, however, applies to any input state with a 
fixed number of photons. Thus, it allows optimization of 
the input state in the presence of an arbitrary amount of 
propagation loss in the two arms of an interferometer. 

With this scheme, the pure state input is now described 
as an (N + 1) x (N + 1) density matrix. After passing 
through the two beam splitters — representing the two 
lossy channels in the two arms of the interferometer — 
the density matrix for the two main modes now consists 
of N + 1 different blocks. Each block represents a given 
number of photons lost. The quantum state of light ends 
up in a mixed state associated with an N + 1 — L dimen- 
sional Hilbert space, if a total of L photons are lost. 

For the arbitrary input state iV^input presented above, 
we find that photon loss can be described by a trans- 
formation to a reduced density matrix whose matrix ele- 
ments are simply given by 
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Here, t, r and t' , r' are transmission and reflection coeffi- 
cients for the fictitious beam splitters in the upper path 
and the lower path, respectively, and C is the number of 
combinations, CJ! = (?). Note that for a given number 
of photons N, the L value labels the block of the reduced 
density matrix, and i, j specify the matrix clement inside 
that block. 

Now we compare a classical interferometer with the 
optimized quantum state approach. Consider a singlc- 
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mode coherent state \a) = e ~ J2T=o ~Jk\\k) as the 
input state to a Mach-Zehnder interferometer (MZI) 
[131 ] . The first beam splitter partitions the state be- 
tween two modes. In the first mode, the state acquires 
a phase shift </> and undergoes a loss of intensity by 
a factor of \t\ 2 . Then, the two beams are redirected 
to the second beam splitter, and photons are detected 
in each output port. The Fisher information, normal- 
ized to the average number of photons n = \a\ 2 , is 
F/n = (4|i| 2 sin 2 60/(1 + \t\ 2 tan 2 6), where 6 describes 
the angle of a y rotation by the first beam splitter. This 
equation provides the Cramer-Rao asymptotic accuracy 
of measurement of the unknown phase shift <fi using a 
classical scheme. We single out two cases. First, we as- 
sume that the first beam splitter has a fixed 50-50 ratio. 
Then in the limit of large loss t — > wc have 



F/f 



2\t\ 



l*l-o 



1+ t 



21*1 



(3) 



Fisher information can be increased by optimizing the 
first beam splitter to compensate for the loss. In this 
case we obtain 



F/n = 
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The optimal quantum input state is now obtained 
numerically. First, a forward problem solver is developed 
using a density matrix approach. An input state is 
written as a density matrix. Phase shifts during photon 
propagation are taken into account by operating with 
e l ^ nA on this density matrix. Next, photon losses are ap- 
plied using Eqs. |T]) and (O, producing a reduced density 
matrix of dimension (N + 2)(N + l)/2. In the last step 
of the forward problem solver, the minimum detectable 
phase sensitivity S<fi is computed from the final density 
matrix. The phase detection is modeled by a final 50-50 
beam splitter followed by two number resolving photode- 
tectors. The joint probability of simultaneously detecting 
mi photons at the first photodetector and 7712 photons 
at the second photodetector is computed as P m — 
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where the label m represents a pair of numbers (mi, 
m-i). Here, XJ\, S is a unitary transformation repre- 
senting a 50-50 beam splitter and /5 ou t is the final 
density matrix obtained after loss. Then, phase sen- 
sitivities are estimated from the Fisher information, 
F, for a single measurement, 86 = 1/y/F, where 

F = p^ m+1)/2 P m (dlnP m /d4>) 2 1. We note 
that, in all of our calculations, we assume a large flux 
of entangled states is to be used, and we normalize our 
results by this flux. 

We optimize the system to find the minimum de- 
tectable phase sensitivity, given fixed losses in the de- 
tection and control arms. For this, a genetic global op- 
timization algorithm is applied to the forward problem 
solver. The parameters to be optimized are the complex 
coefficients Ck\ the optimal sensitivity is necessarily 6- 
independent since a change in 6 can be absorbed into 
the relative phases of c^.. During the numerical com- 
putation of F, we observe that the landscape of F in 
the optimization parameter space possesses several local 
maxima contrary to the convex Fq used by Doner et al. 
i- 

The results of numerical optimization of 86 are pre- 
sented in Fig. O We denote the losses in dB in the de- 
tection and control arms as Ra and Rb, respectively. 
First, in Fig. [U^a) we investigate the overall influence of 
loss in the control arm. One set of simulations is con- 
ducted with equal losses in the detection and control arms 
(Ra = Rb)- We also consider fixed 10 dB loss and fixed 
dB loss in the control arm (R B = 10 dB, R B = dB) 
as loss in the detection arm is varied. In all cases, N = 6 
is assumed. We consistently find that an increase in Rb 
results in higher 5(f>. Thus, one can expect the best phase 
sensitivity to be achieved with the smallest possible loss 
in the control arm. 

In Fig. Ojb) , the numerically optimized phase sensitiv- 
ity is presented as a function of Ra for N = 1, 3, 6, 10, 
and 20. The classical (coherent light) baseline of Eq. © 
is also shown. To compare the quantum results for dif- 
ferent N with each other and with the classical case, we 
rescale the phase sensitivity of the pure quantum state by 
normalizing the Fisher information similarly to Eqs. ([3|) 
and 0: 5cf> = 1/y/F/N = 86\^N. Since 5(f> is obtained 
from the Fisher information per single photon, it is also 
the measure of the synergically-enhanccd phase sensitiv- 
ity per single photon with TV photons acting together. 
For pure quantum states, Fig. [2](b) shows that larger N 
produces smaller 86 for any given amount of loss. The 
iV-dependence of 8<j> is greatest at Ra = dB, and weak- 
est in the limit of extremely high loss, where the lines 
merge together. Coherent light does not show enhance- 
ment with N at any level of loss, and at every level of 
loss, coherent light exhibits worse performance in phase 
sensitivity compared to entangled quantum states. 
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FIG. 2: (Color online) Minimum detectable phase sensitivities 
calculated from the normalized Fisher information: (a) As a 
function of detector arm loss (Ra) with three different losses 
in the control arm (Rb) for N = 6. (b) As a function of Ra for 
Rb = in log scale, (c) As a function of input photon number 
N for Rb = with fixed Ra- Lines represent the results of 
curve fitting using a functional form 1 /N x . In the absence of 
loss, the Heisenberg limit, x = 1, is obtained. For high loss, 
x tends toward 0.5, approaching the shot noise scaling. 



One interesting observation is that 86 in the extremely 
high loss region (Ra >~ 16 dB) becomes AMndcpcndent 
even when using optimally entangled quantum states. In 
other words, the optimal phase sensitivity 86 given by the 
optimal quantum state becomes proportional to N" 1 / 2 
in this high- loss regime, i.e., its scaling with A" in this 
regime is the same as for coherent light governed by the 
shot noise limit. However, despite the same scaling with 
N , the phase sensitivity 86 is still better with entangled 
quantum states than with coherent light. This can be 
explained by the optimal preparation of the initial state. 
As we will see later, the probability amplitudes Ck are dis- 
tributed asymmetrically to generate the smallest possible 
86 for nonzero loss, while coherent light always enters the 
system through 50-50 beam splitter, i.e., it is symmetrical 
between the control and detection arms, Eq. ([3]). With 
coherent light, a similar improvement can be achieved 
by adjusting (or optimizing) the first beam splitter, re- 
sulting in Eq. ([?]). In the latter case, 86 of coherent light 
becomes identical to that of the pure quantum state with 
N = 1. However, we emphasize that, when losses are not 
too high, the phase sensitivity of Af > 1 pure quantum 
states is always better than that of coherent light, even 
with an optimized first beam splitter. 

Figure [2{c) shows the optimal phase sensitivity as a 
function of photon number A" for given fixed Ra values, 
with Rb = 0. Here 86 is shown for each N, with circles 
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FIG. 3: (Color online) (a) Projection of the optimal state on 
NOON, M&M', and GPCS as a function of R A , for N = 4. (b) 
The optimal input state composition. The vertical axis shows 
coefficients \ck\ 2 of the optimally entangled input state. Both 
figures clearly demonstrate that the optimal state changes 
from NOON-type to GPCS-type as loss increases, with the 
crossover occurring at approximately 5 dB loss. 



indicating N = 1, 3, 4, 6, 10, and 20. The lines are 
drawn by curve-fitting to a power law, 8<p ~ 1/Y X . From 
this, we find that 8<j) is well represented by 1/N, 1/Y - 68 , 
l/Y - 517 , and 1/iV 5 for R A = dB, 2.88 dB, 10.02 
dB, and 21.86 dB, respectively. This result provides an 
overall view of how the phase sensitivity changes from the 
Heisenberg limit to the shot noise limit with increased 
loss, i.e., it follows a power law with \/N x where 1/2 ^ 
x s$ 1. 

To characterize the optimal state, we use three classes 
of well-defined states: NOON, M&M', and a two- 
mode SU(2) coherent state, often called the General- 
ized Perelomov Coherent State (GPCS). GPCS is defined 



as \GPCS) 
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where aj and a\ are creation operators in the two modes, 
and a and (3 are two real parameters 15(. In particular, 
M&M' is the first class of path entangled states shown 
analytically to have robustness to photon loss. It is in- 
teresting to see how the true optimal state may differ 
from the M&M' state in lossy environments. Charac- 
teristics of the optimal state are presented in Fig. [3] for 
a fixed Rb = dB. The similarity of the optimal states 
with each of these three benchmark states is measured 
by the squared overlap between the optimal state and 
the benchmark state. The results show that the optimal 
state is closest to the NOON state for low loss. As loss in- 
creases, the NOON state portion gradually decreases and 
the optimal state becomes closer to GPCS than to NOON 
at around 5 dB of loss. The degree of similarity between 
the optimal state and the M&M' state is rather low for 
loss smaller than 6 dB. For every value of loss, GPCS 
is higher than M&M'. Figure [^b) shows how the input 
amplitudes of the optimal state are arranged for different 
loss levels. In the lossless case, we have the NOON state. 
As loss increases, the optimal state is reshuffled and ac- 
quires an asymmetric shape. This serves as critical in- 



formation for achieving a highly sensitive intcrfcromctric 
system. Based on the results shown, it is obvious that 
generating such optimal input states should be the first 
consideration in the development of an intcrferometric 
sensor using entangled photons. 

In summary, we have performed unconstrained opti- 
mization of a lossy two-mode interferometer. Wc con- 
clude that input NOON states are optimal for nearly zero 
loss [1], and that finite-photon number two mode coher- 
ent states are optimal — with shot-noise sensitivity — 
for large loss. Our results suggest that, if sensitivity is 
the only metric of success, ordinary coherent input state 
interferometry is best for high loss. This leaves open 
super-sensitive schemes employing squeezed light at the 
detector [lH ] or super-resolving schemes employing pho- 
ton number resolving detectors [13] ■ 
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